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ABSTRACT 

Peculiar velocities arise from gravitational instability, and thus are linked to the 
surrounding distribution of matter. In order to understand the motion of the 
Local Group with respect to the Cosmic Microwave Background, a deep all-sky 
map of the galaxy distribution is required. Here we present a new redshift com- 
pilation of 69 160 galaxies, dubbed 2M++, to map large-scale structures of the 
Local Universe over nearly the whole sky, and reaching depths of K < 12.5, or 
200h~^ Mpc. The target catalogue is based on the Two-Micron- All- Sky Extended 
Source Catalog (2MASS-XSC). The primary sources of redshifts are the 2MASS 
Redshift Survey, the 6dF galaxy redshift survey and the Sloan Digital Sky Sur- 
vey (DR7). We assess redshift completeness in each region and compute the 
weights required to correct for redshift incompleteness and apparent magnitude 
limits, and discuss corrections for incompleteness in the Zone of Avoidance. We 
present the density field for this survey, and discuss the importance of large-scale 
structures such as the Shapley Concentration. 



1 INTRODUCTION 



Peculiar velocities remain the only method to map 
the distribution of dark matter on very large 
scales in the low redshift Universe. Recently, sev- 



revious full- sky gala xy redshift catalogues like PSCz 
fl2000l) or the Two-Micron- All-Sk> 



(|Saunders et al. ,,-. 
shift Survey (IHuc 



IJ or tne i wq-lVlicron-Aii-bky Red- 
. et al. 1 12OO5I : lErdogdu et all l2006al : 



Watkins et al.ll2009l:l] 


^javaux et al...201Q,;„Kashlinskv et al.l 


2OIOI: iFeldman et al. 


I2OIOI) of the mean or "bulk" flow 



^ Mpc suggest a high veloc- 
Mpc volume with respect 



on scales larger than lOO/i 
ity of our local ^ 100 
to the Cosmic Microwave Background (CMB) frame. 
In the standard cosmological framework, peculiar ve- 
locities are proportional to peculiar acceleration and so 
one expects the bulk flow to arise from fluctuations in 
the distribution of matter, and hence presumably of 
galaxies, on very large scales. Another statistic for mea- 
suring such large-scale fluctuations is the convergence 
of the gravity dipole as a function of distance. How- 
ever, the rat e of convergence has been a subject of re- 



iHuchra et al.ll2011 . 2MRS). We present here a new cata- 
logue called the 2M++ galaxy redshift compilation. The 
photometry for this compilation is based is on the Two- 
Micron- All-Sky-Survey (2MASS) Extended Source Cata- 
log dSkrutskie et al. 200fi, 2MASS-XSC). We gather the 
hi^h-quality redshif t s from the 2MASS redshift survey 
(iHuchra et al.ll2005l : lErdogdu et alll2006al : iHuchra et all 
l201lh limited to K = 11.5, the 6dF galaxy redshi ft survey 
Data Release Three (GdFGRS IJones et aLl lioOol) and the 
Sloan Diffltal Sky Surv ey Data Release Seven (SDSS-DR7 



lAbazaiian et al.ll2009l ) 



cent d e bate (iKocevski Ebelind l2006a 
l2006al lbl: lLavaux et al.ll2010l : lBilicki et al 



Erdogdu et al.l 



I2OIIL and refer- 



ences therein) . A closely-related topic is the gravitational 
influence of the Shapley Concentration, the largest con- 
centration of galaxy clusters in the nearby Universe. It 
is therefore important to have catalogues that are as full 
sky and as deep as possible to understand whether the 
distribution of matter in the nearby Universe may explain 
the above-mentioned results. 

It is already possible with currently available data 
to build a redshift catalogue significantly deeper than 



A summary of this paper is as follows. In Section [J] 
we describe the steps in constructing the 2M++ red- 
shift galaxy catalogue: source selection, magnitude cor- 
rections, redshift incompleteness estimation and correc- 
tion, the luminosity function (LP) estimation and the 
final weight computation. In Section [S] we discuss the 
Zone of Avoidance (ZoA) in our catalogue, and how its 
effects can be mitigated. In Section |4l we define groups 
of galaxies and check some of their overall properties. In 
Section [5l we compute and analyse the density field, pre- 
senting maps of the Supergalactic plane and three cluster 
density and velocity profiles. Section [6] summarizes our 
key results. 
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2 CATALOGUE CONSTRUCTION 

In this Section, we describe the construction of the 2M++ 
galaxy redshift catalogue from the different data sources. 
First, in Section 12.11 we describe the source data sets 
that form the basis of our catalogue, as well as the pri- 
mary steps in the construction of the 2M++ catalogue. 
We then present the methodology used for merging these 
different sources. In Section \T2l we describe the correc- 
tions applied to apparent magnitudes to homogenize the 
target selection. In Section [231 we test and apply the red- 
shift cloning procedure to our data to increase the overall 
redshift completeness. We then estimate redshift com- 
pleteness (Section I2.4.2P and present the number counts 
of galaxies as a function of redshift (Section [23|). Finally, 
we compute the LF of our sample in Section 12.61 and 
compute the total weights to apply to each galaxy in 
Section [221 



2.1 Source datasets and construction procedure 

Our ca talogue is based on the Two-Micron- All-Sky- 
Survey fekrutskie et al.ll2006l . 2MASS) photometric cat- 
alogue for target selection, w hich has very hig h com- 
pleteness up to Ks = 13.2 (jCole et al.1 l200lh . Here- 
after, for brevity, we use K in place of Ks. As noted 
above, we will be using redshifts from the SDSS- 
DR7, the 6dfGRS and the 2MRS. In addition to 
these main sources, we have gathered additional red- 
shifts from a number of other so u rces (jSchneider et al.l 
'1990"; de Vaucouleu rs et al.l 
Huchra et al. 1995: 



T991I : 'Binggeli et al.1 Il993l: 
Falco et al. 

I [l999;" Conselice et al.l 



l200li r iRines et alTl2003l : iKoribalski et al.ri2004l ) through 
NED queries Q Due to the inhomogeneity of the tar- 
get selection between the different redshift surveys, we 
think that it is more appropriate to define a new tar- 
get selection rather than using existing target databases 
from the above surv eys. We used of the NYU-VAGC 
(|Blanton et al.l I2OO5I ) catalogue for matching the SDSS 
data to the 2MASS Extended Source Catalog (2MASS- 
XSC). The NYU-VAG C provides the S DSS survey mask 
in MANGLE format ([Hamilton Tegm ark 20M)Q We 
sampled the mask on a HEALPix grid at A^side = 512 
('^10 arcminutes resolution). This angular resolution cor- 
responds to ^lh~^ Mpc at ^300/i~^ Mpc. Because ulti- 
mately we will be smoothing the density field on scales 
of ^4h~^ Mpc, the mask has sufficient resolution for our 
purposes. Additionally, we filter out from our target se- 
lection the extended sources that are known not to be 
galaxies 

We aim to limit 2M++ at K ~ 12.5 regions of the 
sky covered by SDSS or by 6dF. The exact cut depends 
on the adopted definition for the magnitude. As we want 



1 The NASA/IPAC Extragalactic Database (NED) is oper- 
ated by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aeronautics and 
Space Administration. 

^ We use the file named lss_combmask.dr72.ply, which gives 
the geometry of the DR72 sample in terms of target selection 
with bright stars excised. 

^ We require that the VlSUAL_CODE is not equal to two. 



to retain as much as possible information from the shal- 
lower 2MRS catalogue, we opt to follow closely the mag- 
nitude used by 2MRS for target selection. We define as 
K2M-\-+ the magnitude of a galaxy measured in the Ks 
band, within the circular isophote at 20 magarcsec"^, af- 
ter various corrections as described below (Section 12. 2p . 
Several of the steps taken to build the catalogue are de- 
scribed in greater detail in the following Sections. We now 
outline these steps: 

(i) We import the redshift information for 2MASS- 
XSC galaxies from the NYU-VAGC for SDSS-DR7, the 
6dF-DR3, and from the 2MASS Redshift Survey. 

(ii) We correct for small-scale redshift incompleteness 
(arising from fibre collisions) by 'cloning' the redshifts of 
nearby galaxies (Section 12. 3p . 

(iii) We correct the apparent magnitudes for Galactic 
dust extinction (Section l2.2p . 

(iv) We use the redshift to correct for galaxy evolution- 
ary effects and aperture corrections fSection l2.2p . We call 
this magnitude i^2M++. At those magnitudes, we assume 
that the photometric completeness is one at Galactic lat- 
itudes higher than 10°. 

(v) We compute two sets of galaxy samples: a tar- 
get sample with K2M++ < 11-5 in regions not covered 
by 6dFGRS or SDSS, and a target sample limited to 
K2M-\-\- < 12.5 in regions covered by SDSS or 6dFGRS. 

(vi) We estimate the redshift completeness as a func- 
tion of position on the sky for in each of these regions. 

(vii) We place galaxies in groups and clusters using a 
percolation algorithm. 

(viii) We compute the Schechter parameters of the LF 
of the combined catalogue (Section |2.6|) . We use this LF 
to compute the weights to apply to each of the observed 
galaxies to take into account the unobserved ones (Sec- 
tion [23|. 

In future work we will update the estimated distances 
for the galaxies using reconstructed velocity field and re- 
execute step (viii) to update the corrections. The detail of 
this procedure will be discussed in a later paper (Lavaux 
& Hudson, 2011, in preparation). 



2.2 Apparent magnitude corrections 

We describe in this Section the corrections that must be 
applied to apparent magnitudes to mitigate the effects 
of cosmological surface brightness dimming. Galactic ex- 
tinction and stellar evolution. We choose to use a defini- 
tion of the magnitudes for target selection that is related 
to the one used for defining magnitudes in the 2MRS. 
This ensures that the final completeness is maximized in 
in the parts of the sky where only redshifts from 2MRS 
are available. 

The absolute magnitude M at redshift zero of a 
galaxy may be written as 

M = m- Ak{1, b) - k{z) + e{z) - Dl{z) (1) 

with m the apparent magnitude, Ak{1^ b) the absorption 
by Milky Way's dust in the direction (/,6), k{z) the k- 
correction due to the redshifting of the spectrum, e{z) is 
the correction for evolution of the stellar population and 
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Figure 1. Error distribution due to the redshift cloning procedure - We give here the computed error of the redshifts of either 
SDSS or 6dF redshift catalogue. We removed the redshift information of half of the objects in these catalogues and tried to recover 
them using the cloning procedure. The difference is plotted as an histogram in the two plots above. The overlaid continuous curve 
correspond to a Cauchy-Lorentz distribution with a width equal to 2.7 Mpc. 



Dl is the luminosity distance. We convert the redshifts 
into luminosity distances assuming a ACDM cosmology 
with a mean total matter density parameter Qm = 0.30, 
and a Dark Energy density parameter Qa = 0.70. All ab- 
solute magnitudes are computed assuming Ho = 100 km 
s~^ Mpc~^ 

The absorption in Ks band is related to the extinc- 
tion E b-v estimated using the maps of ISchlegel et al.] 
(|l998l ) by the relation 

= 0.35^(s-F)(/,&) (2) 

where the constant of proportionality is obtained from 
the relation between absorption in K band and in V band 
(|Cardelli et al.lll989l ). 

The adopted /c-correction is 

k{z) = -2.lz (3) 

from iBell et al.1 (|2003al ). Finally, the evolutionary correc- 
tion is 

e{z) = 0.8z (4) 

also from lBell et~all (l2003al ). 

The magnitudes adopted in this work are circular 
aperture magnitudes defined within a limiting surface 
brightness of 20 mag per square arcsec. However, vari- 
ous effects will cause not only the observed magnitude 
to change, but also the observed surface brightness. As 
the surface brightness of the galaxy profile drops, the 
isophotal aperture will move inwards and so the aperture 
magnitude will drop. 

The surface brightness will depend on redshift via 
the usual (1 + cosmological dimming effect as well 
as the extinction, /c— correction and evolutionary effects 
described above. Therefore the correction is 

ASB = SB(z = 0) - SB(z) 

= -10 log(l + z)- h) - k{z) + e{z) . (5) 

Note that the k-\- e corrections have opposite sign to the 



cosmological surface brightness dimming, and so there is 
some cancellation of these effects. However, the cosmo- 
logical term still dominates, so the net effect is that as 
galaxies are moved to higher redshift their surface bright- 
ness is dimmer. 

By simulating simple Sersic profiles, we have esti- 
mated how much the aperture magnitude changes as a 
result of surface brightness dimming. For a typical 2M++ 
galaxy with Sersic index n = 1.5 and mean surface bright- 
ness within the effective radius of (/Xe) = 17.5, we find 
that the correction to the magnitude due only to a change 
in aperture radius can be approximated by 0.16 ASB 
where ASB is the correction in the surface brightness. 
This term is only the shift in magnitude due to the shift 
in isophotal radius, and does not include the "direct" ef- 
fect on the magnitude itself due to extinction and k -\- e 
corrections. Thus the total effect is: 

K2M++ = K2o,c + 1.16[-Ak(/, b) - k{z) + e{z)] 

-0.16[101og(l + z)] (6) 

Note that this is close, but not identical, to the 2MRS 
corrected magnitude. 

In some cases, only the magnitude K2o,e, derived 
from adjusting an ellipsoidal Sersic profile, is available. In 
those cases, we have computed the corresponding K2o,c 
using the following relation, obtained by fitting on the 
galaxies for which the two magnitudes were available: 

K2o,c = (0.9774 ± 0.0005)ii^2o,e + (0.288 ± 0.006). (7) 

The residual of the fit has a standard deviation equal to 
0.11. We also use this relation whenever the predicted 
K2o,c and the actual K2o,c from 2MASS-XSC differs by 
0.22 and calculate the K2o,c from K2o,e- We have used 
this relation for 7% of galaxies, both in the target and 
the final redshift compilation. 
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2.3 Redshift cloning 

Within the 6dF and SDSS regions, there is small-scale in- 
completeness due primarily to fibre collisions. To improve 
the redshift coverage of the catalogues we "clone" red- 
shifts of nearby galaxies within each survey region. This 
procedure, wh i ch is related to another one described in 
iBlanton et all (|2005l ). is as follows. Consider two targets 
Ta and T^. If Ta does not have a measured redshift and 
Th has one, and furthermore Th is the nearest target of 
Ta with a angular distance less than e, we copy the red- 
shift of Th to Ta. e is determined by the angular distance 
between two fibr es of the measurin g; instrument, which is 
e = 5^7 for 6dF (I Jones et"alll2004l ) and e = 55'' for SDSS 
(jBlanton et al.ll2003h . We refer to redshifts cloned in this 
way as "fibre-clones". 

To assess the errors on redshifts for the fibre-clones, 
we randomly split the set of galaxies which have a mea- 
sured redshifts in two sets <Skeep and <Stest- We mark the 
galaxies belonging to <Stest as having no redshift. We then 
apply the fibre cloning procedure to these galaxies. 

The result of this test is shown in Figure [1] 
for both SDSS galaxies and 6dF galaxies. We note 
that the Cauchy-Lorentz distribution with width W — 
2.7h~^ Mpc is a good fit to the central part of the two 
distributions. We used the formula 



P(e) 



nW l + {e/Wy 



(8) 



for the modelled probability distribution function in both 
panels. We checked that a Gaussian distribution man- 
ages only to fit the central part of the distribution and 
is less adequate than a Cauchy-Lorentz distribution. The 
fibre-clones are given a redshift error of 9 X 10"^, which 
corresponds to ^ 2.7h~^ Mpc at redshift z = 0. 



2.4 Redshift survey masks and completeness 

Because of the different redshift catalogues used in 
2M++, we will separate the full sky into the follow- 
ing regions: K2M++ < 11.5 (2MRS); K2M++ < 12.5 
(6dF/SDSS) or regions with insufficient redshift data. We 
begin by assigning preliminary "masks", then measure 
the redshift completeness and then assign final masks 
based on a completeness limit of 50%. We now describe 
these steps in more detail. 



2.4- i Preliminary mask selection 

The prelim i nary mask for the 6dF is as given in 
I Jones et al.1 (120041 ) : \b\ > 10° and in the southern hemi- 
sphere ^ < 0°. 

For SDSS, from the full DR7, we select only the most 
homogeneous and contiguous portion of the redshift sur- 
vey. To obtain the geometry corresponding to this por- 
tion, we use the mask computed numerically from the 
MANGLE file, and impose an additional constraint on 
the positions: we keep only galaxies within the region 
90° < a < 250°, to which we add another region at 
250° < a < 270° and ^ < 50°. This selection retains 
the major contiguous piece of the SDSS in the North- 
ern Galactic cap, while removing the Southern Galactic 




-90 

Figure 2. Molleweide projection of the SDSS spectroscopic 
mask in Equatorial coordinates. We removed unconnected 
parts from the original mask. Note the presence of small holes 
in the mask due to both the presence of stars and not exact 
reconnection of the SDSS plates. This corresponds to an inter- 
section of the geometry described by lss_geometry.ply and 
our selection criterion described in Section [274] 




2MRS 



SDSS 



6dFGRS 



Figure 3. A Molleweide projection in galactic coordinates 
showing the preliminary masks corresponding to the differ- 
ent redshift catalogues: the 2MASS Redshift Survey (red); the 
6dF galaxy redshift survey (green), and the Sloan Digital Sky 
Survey DR7 (blue). Regions with no redshift data are shown 
in grey. 



strips and the small disjoint piece with a ^ 260. After 
these cuts, we still retain 90% of the area covered by the 
complete SDSS-DR7 "Legacy" spectroscopic survey. 

The preliminary SDSS mask is shown in Fig. [2] This 
mask has a relatively simple geometry and is contiguous, 
except for the presence of very small holes that are due 
either to stars or to imperfect overlap of the SDSS plates. 
These imperfections represent about 8% of the SDSS sur- 
face area. 

For 2MRS, the initial mask is the whole sky with 
\b\ > 5°, except in the region —3 0° < / < +30°, wher e 
the Galactic latitude is |6| > 10° (jErdogdu et al.ll2006al ). 
and excluding the regions covered by 6dF or SDSS. We 
refer to this region as 2Mx6S. 

The combination of the three survey masks is shown 
in Fig. [3l in Galactic coordinates. The grey area near the 
galactic plane is covered by none of the surveys. There 
is a small overlap between the SDSS and the GdFGRS in 
the north galactic cap (in cyan). 
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2.4-2 Redshift completeness 

In order to determine LFs and weights needed for the den- 
sity field, it is first necessary to assess the completeness 
of the redshift catalogues on the sky. We will estimate 
the redshift completeness in some direction of the sky 
for two different magnitudes cuts: X2M++ < 11-5 (for 
the 2Mx6S region) and K2M++ < 12.5 (for 6dF/SDSS 
regions). The maps of these completeness are given in 
Fig. m The completeness of the 2MRS is quite homoge- 
neous and only drops close to the Galactic plane. The 
6dF survey is mostly homogeneous except at the loca- 
tion of the Magellanic Clouds. The SDSS completeness is 
quite homogeneous, and remains at a level of about 80% 
in the whole contiguous region. 



2.4-3 Final masks and the 2M++ catalogue 

Based on the above analysis, we reject from 6dF and 
SDSS those regions where the completeness at < 
12.5 is less than 50%, and assign these areas to the 2Mx6S 
and limit the magnitude there to 11.5. As a result, we 
note that while the 6dF and SDSS masks are contigu- 
ous, the 2Mx6S mask is not. We show in the lower left 
panel of Fig. O the mask corresponding to the footprint 
of the 2MH-H- compilation. The black corresponds to re- 
gions where the survey has a completeness higher than 
50% at the limiting magnitude of K2M++ — 12.5 (the fi- 
nal 6dF and SDSS masks). The grey area represents the 
same but for a limiting magnitude of = 11.5 (the 

final 2Mx6S mask). In white, we show the parts of the 
sky where either there is no redshift information or target 
galaxies were not present. 

The Zone of Avoidance is clearly visible. There are 
also important unobserved patches at K2M++ < 12.5 in 
the southern Galactic hemisphere at the locations of the 
Magellanic Clouds. The other white patches in the south- 
ern hemisphere are mostly related to local higher absorp- 
tion by dust in the Milky Way. 

Thus, in summary, the final 2M++ catalogue is de- 
fined as all galaxies in 2Mx6S with K2M++ < 11.5, or in 
6dF or SDSS with K2M++ < 12.5, and contains 69 160 
galaxies with redshifts (including fibre-clones). Table [T] 
summarizes the statistics and completeness for the differ- 
ent regions of the 2M++ compilation. We note that the 
2M++ compilation redshifts are nearly 90% complete, 
and so redshift completeness corrections are small. 



_c5 
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Figure 5. Redshift distribution of 2M++ galaxies. We show, 
in grey histogram, the number of galaxies within each redshift 
bin 5z = 0.00150. The cumulative number of galaxies at a red- 
shift less or equal to z is given by the red curve. The predicted 
number of galaxies given by our fiducial LF given at the end 
of Section I2. 61 is shown in solid green for the cumulative num- 
ber and in solid blue for the number of galaxies in each bin of 
the grey histogram. The LF has been fit using a subset of the 
catalogue for which 5,000 km s~-'^< cz <20,000 km s"-*^ and 
-25 < M < -21. 



Distribution for K2M++ < 11-5 
Distribution for K2M++ < 12.5 




0.06 0.08 

z (redshift) 

Figure 6. . Redshift galaxy distribution in 2M++ for 
K2M++ < 11-5 (filled histogram) and K2M++ < 12.5 (un- 
filled histogram). 



2.5 Redshift number density of galaxies 

Within the three regions outlined above, there are a total 
of 69 160 galaxy redshifts (including fibre-clones). Fig. [5] 
shows a histogram of all redshifts, as well as the cumula- 
tive counts starting from redshift z = 0. Conservatively, 
the catalogue appears totally complete up to 2: = 0.02 
(- 60h~^ Mpc). This is due to our use of the 2MRS for 
one part of the sky. 

In Fig. [6l we compare quantitatively the counts in 
the 2MRS region with K2M++ < 11.5 with the 6dF/SDSS 
regions with K2m-\-\- < 12.5. Because our magnitude cor- 
rections are not precisely equivalent to those used for the 
2MRS catalogue, the increase of the magnitude cut to 



K2M++ = 12.5 is not strictly equivalent to using only the 
6dF and SDSS spectroscopic data but also includes a few 
2MRS galaxies. Nonetheless, the increase in the magni- 
tude cut correspond mostly to the the sky portions cov- 
ered by both the SDSS and the 6dF. The deeper redshift 
data allows us to better probe more distant large-scale 
structures, particularly in the redshift 0.02 ^ ^ ^ 0.05. 
For example, the feature in the redshift distribution at 
z ^ 0.05 corresponding to the Shapley Concentration is 
not present for the subcatalogue i^2M++ < 11.5, while it 
is clearly seen in the 6dF subcatalogue {K2m+-\- < 12.5). 
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Figure 4. The 2M++ compilation in Galactic coordinates. The top panels give the redshift completeness of 2M++ for the limiting 
magnitude K2M++ = 11.5 (left panel) and K2M++ = 12.5 (right panel). The grey area corresponds to regions with no targets and 
no redshifts. In the bottom left panel, we show the final combined mask. The regions in white have been to completely excluded 
because of very low completeness. The regions in grey (black respectively) have a redshift completeness higher than 50% at a 
limiting magnitude K2M-\--\- ^ 11-5 (i^2M++ ^ 12.5 respectively). In the bottom right panel, we show all galaxies with redshift 
included in the 2M++ compilation. Each galaxy is color-coded according to its redshift distance, blue for the nearest and red for 
the farthest. 



Table 1. Summary statistics for the primary regions in the 2MH — h compilation 



Region Area Nm<mii^ Nz c 

2Mx6S 11.5 13 069 9 419 9 016 0.96 

6dF 12.5 17 041 46 734 42 442 0.91 

SDSS 12.5 6 970 20 333 17 702 0.87 

None - 4 172 - - - 

Total - 37 080 76 451 69 160 0.90 



Note that the regions are counted exclusively. We have not enforced the sample to have a local redshift completeness higher than 
50%, resulting in a total number of redshifts higher than in the final catalogue. 



2.6 Luminosity function 

2.6.1 Method 

In order to correct for selection effects due to magni- 
tude limits, it is first necessary to measure the LF. We 
take into account the redshift completeness to measure 
the LF of galaxies in the combined catalogue. We use 
a m odified version o f the likelihood formalism used to 
find Schechter (1976) function parameters, as described 
by bandage et al. (197^). We assume that evolutionary 
effects on the luminosity of galaxies have been accounted 
for by Eq. (|6|). The parametrization adopted is the usual 



Schechter function: 




with n* the density normalization constant, L* the char- 
acteristic luminosity break, or equivalent ly in terms of 
absolute magnitudes 

$(M) = 

0.41og(10)n*10°-^(^+-)^^*-^) exp (_ioO-4(^*-^)) 

= n*$o(M), (10) 
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Figure 7. Galaxy LF estimates. The left-hand panel shows the non-parametric LF estimated using the 1/Vmax method is shown 
by the solid lines for the regions covered either by 2MRS, 6dF, SDSS or all together. For these plots we use data from 750 km s"-*^ 
to 20 000km s"-*^. The dashed line shows the parametric LF using the likelihood method of Section [2?6] for galaxies with absolute 
magnitudes in the range —25 < M < —21, for redshift distances 5 000 km s"-*^ to 20 000 km s"-*^. The error bars reflect only the 
uncertainties in galaxy counts and do not include cosmic variance effects. The right-hand panel shows the difference between the 
1/Vmax LFs and the fltted parametric LF. 



with M* the characteristic absolute magnitude break 
in the Schechter function. Above, we have introduced 
$o(M), which is the unnormalized Schechter function. 
We model the probability of observing a galaxy of abso- 
lute magnitude Mi given its redshift Zi as 



P(M^\z^,a,M\n%c) 



c{Mi,Ui,di)^o{Mi) 
/^-^^ c(M, u,,d,)MM) dM 



(11) 



with Mmin, Mmax the maximum absolute magnitude 
range from which the galaxies were selected in the cat- 
alogue, c{M,Ui,di) the completeness in the direction Ui 
of the object i, at the absolute magnitude M, di the lu- 
minosity distance of the galaxy i at redshift Zi. This ex- 
pression is simplified using our assumption that redshift 
incompleteness c(M, u, r) may be modelled by two maps 
at two apparent magnitude cuts. c(M, r) is thus 



c{M,u,r) 



ct{u) if M + 51ogio (to^) < ^6 

Cf{u) if rrib < M + 5 log^o (to^) < 




otherwise, 

(12) 

with rrih = 11.5 and irif = 12.5. The expression of the 
probability (|11|) may thus be newly expressed as 



P{M^\z^,a,M\n,c) = 



c{M^,Ui,d^)^o{M^] 
/(di,'Zii,Mmin,Mma> 



(13) 



with 



f{r,f, Mmin, Mm..) = ct{r)rl 



mm ) ^^-^ 

(14) 

the normalization coefficient for the direction f at dis- 
tance r, and rio defined as the distance in units of 10 pc. 



fmb,rio) + 



In the above, we have also used the function F^ '°^,Mmax 
defined as 



'."^ ^ (m,rio) = 

r( -i I -1 /-v0.4( M* — min(max( 

inc ( 1 + a, 10 ^ 



M(m,rio),M^in),Mmax))^ 



(15) 

with the absolute magnitude 

M(m,rio) = m - 51ogio(rio) (16) 
and rinc(tt,2/) the incomplete Gamma function 



I a—1 —X 

X e 

y 



dx . 



(17) 



We write the total probability of observing the galax- 
ies with intrinsic magnitude {Mi} and redshift {zi} given 
the Schechter LF parameters: 

PaMi}\{ci},{zi},(^,M*,n'-) = 

N 1 

-'^galaxies 

Yl P{M^\z^,a,M\n\ci). (18) 

Using Bayes theorem, we now estimate the most likely 
value taken by a,M* assuming a flat prior on these pa- 
rameters. 

The normalization constant n* is determined us- 
i n^ th e minimum variance estimator of IPavis Hucfiral 
(|l982l ), but neglecting the effects of cosmic variance on 
the weights by setting J3 = 0. While our estimate may 
be biased relative to the galaxy mean density outside the 
catalog, it is less noisy than the optimal case. The esti- 
mate also corresponds better to the density in the piece 
of Universe that we consider than the density correspond- 
ing to the optimal weighing. Our choice leads also to a 
simplification of the mean density as the total number of 
galaxies divided by the effective volume of 2MH-H-. Con- 
sequently, for a survey limited in the absolute magnitude 
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range [Mmin, Mmax] and with volume V we compute the 
mean density of galaxy n by the following equation 



galaxies 



/^d2r/(r,f,Mmin,M^ 



with a standard deviation only from Poisson noise 

axies 
^ -^^ffal axies 



(19) 



(20) 



because we have set Js = 0. We then convert n into n* 
using 



(21) 



Im^.^^ MM) dM' 



Similarly it is possible to define the luminosity density 

L = n*lO°-^^^0-'^*^r(2 + a) X (IL©), (22) 

with rf Qj) — i inc (a, 0). The luminosity density is less sen- 
sitive than the number density to fluctuations in a. 

To determine the LF parameters, we select a subset 
of the galaxies in our catalogue. We have defined the 
subset by the joint conditions: 

- Galaxies must have a redshift z such that 5,000 km 
s~^< cz < 20,000 km s~^. The lower limit reduces the 
impact of peculiar velocities on absolute magnitude esti- 
mation, which is derived using redshift s in the CMB rest 
frame. By limiting to cz < 20, 000 km s~^, we avoid more 
distant volumes with high incompleteness. 

- The absolute magnitude estimated from the red- 
shift in CMB rest frame is within the range [Mmin = 
— 25, Mmax = —21]. As mentioned later in this Section, 
this magnitude selection removes the bright objects that 
do not seeni to fo llow a Schechter LF (as also discussed 
I Jones et al.1 120061 ). 

Absolute magnitudes are determined with H = lOO/i km 
s~^ Mpc~^with h = 1 and we have assumed a flat 
ACDMcosmology = 0.30 and Qa = 0.70. We do 
not distinguish between early-type and late-type galax- 
ies, and so fit both populations with a single parameter. 



2.6.2 Results 

The derived LF parameters are summarized in Table [2]for 
our default choice of cuts discussed above as well as for 
other choices that we discuss below. The error-bars are 
given at 68% confidence limit, estimated using a Monte- 
Carlo-Markov-Chain method. 

Fig. [7] shows the LF for our default cuts in the CMB 
rest frame. We also show, for the entire data set and for 
each sub catalogue, the non-parame tric LFs estimated us- 
ing t he unbiased 1/Vmax method ()Schmidtlll968l : iFeltenI 
1 19761 ). Note that for the 1/Vmax LFs the volume and mag- 
nitude limits are different than for the parametric fit, 
which explains that the fitted parametric faint-end slope 
is not a good fit to the 1/Knax- in the range [—21, —19]. 
In the left panel, we give the LFs and in the right panel 
the ratio between the estimated LFs and the best fit. 



2. 6. 3 Discussion and comparison with previous results 

Table [2] also lists LF parameters from previous 2MASS 
studies. Our fitted LF parameters are in agreement 
with previous studies of the K-band LF (|Bell et al.l 
l2003bl : lEke et al.ll2005l) but are somewhat different tha n 
those found bvlKochanek et al.1 diooH) . ICole et al.1 (|200lh . 
iHuchra et alT(l2005l ) and I Jones et all (|2006l ) . 

The derived LF parameters are sensitive to a number 
of systematic effects: the magnitude range used, the rest- 
frame used for the redshifts, and the fitting method itself. 

The Schechter function itself appears not to perfect 
fit over the whole range of magnitudes. Consequently, 
the fitted parameters depend in the magnitude (and dis- 
tance) range of the galaxies used in the fit. Our de- 
fault minimum distance r > 5000 km s~^ corresponds 
to Mk < -21 for = 12.5. However, the 1/Knax 

method seems to indicate an infle ction in the LF at M k ^ 
-21. This bend is also seen bv ISell et al.1 (l2003bl ) and 
lEke et al.1 (l2005l ). Indeed lBell et al.1 ('2003b) attempted to 
fit the part at Mk > — 21 with a po wer-law instead of a 
Schechter functi on.) Several studies (|Biviano et al.lll995l : 
lYagi et aL I I2OO2I ) have noted a dip in the LF of cluster 
galaxies at a similar location (approximately 2 magni- 
tudes below M*), although othe r studies su ggest that it 
is a flattening rather than a dip (|Trentham|[l998t ) . In any 
case, it seems clear that the choice of mag nitude range 
will a ffect the Schechter LF parameters. In Ijones et al.l 
(|2006l ) and lCole et aD (I2OOII ). the magnitude range used 
in the fit is fainter than our default. 

A second issue, which arises when using galaxies with 
very low redshifts, is the choice of flow model or rest- 
frame redshifts. Very nearby galaxies are likely to share 
the peculiar velocity of the Local Group (LG), so the red- 
shift in the LG frame is a better proxy distance than the 
CMB-frame redshift. For better understanding of the de- 
pendence of our results on both local flows and clustering, 
we have fit the parameters of the Schechter function in 
two rest frames (CMB or LG). We find that, for samples 
extending to Mk —17, the faint-end slope a steepens, 
but only by 0.06. 

Finally, the magnitudes, correct and the fitting 
method itself are probably the most important system- 
atics. 



(i) We note t hat t he studies of ICole et al.1 (I2OOII )- 
lEke et al.1 (l2005l ) and iBell et al.1 (l2003bl) are~ based on 



Kron magnitudes, and that of ljones et al.l (|2006l ) is based 
on total magnitudes, leading to a possible difference i n 
M* of 0.20=b0.04, as discussed bv lKochanek et al.1 (I2OOII). 

(ii) Another notable difference is that Ijones et al.l 
(|2006l ) have tried to integrate the effect of uncertainties 
on the determination of magnitudes, which we do not do 
here. 

(iii) iBell et al.1 (l2003bl ) matches SDSS redshi fts to both 
the 2M ASS XSC and the PSC catalogues. iBell et al.1 
(|2003bl ) argues that selection effects bias the raw 2MASS 
LF compared to the true LF. However, whereas those 
authors were interested in, for example, the total stellar 
mass density in the nearby Universe, our goal is rather 
a consistent magnitude system coupled with uniform se- 
lection across the sky. Since our primary method will be 
to weight by luminosity, the small missing contribution 
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Table 2. Summary of K-band Schechter LF parameters from this paper and selected results from the literature. Magnitude ranges 
with a ~ are estimated, is in units of {10~^h^ Mpc~^) and the luminosity density L is in units of lO^h LqMpc"^, assuming 
Mk,q = 3.29 





r L cilllt: 


Magnitude 








71* 




Kochanek et al. (2001) 




[-26 


-20] 


[2000; 14000] 


-1.09 ±0.06 


-23.39 ±0.05 


1.16 ±0.10 




Cole et al. (2001) 






26; -20] 


? 


-0.96 ±0.05 


-23.44 ±0.03 


1.08 ±0.16 




Bell et al. f2003b) 






25; -18] 


? 


-0.77 ±0.04 


-23.29 ±0.05 


1.43 ±0.07 




Eke et al. ^2005) 






25; -20] 


? 


-0.81 ±0.07 


-23.43 ±0.04 


1.43 ±0.08 




Huchra et al. (2005) 






28.5; -16] 


? 


-1.02 


-23.4 


1.08 




Jones et al. (2006) 




[-28.85; -15.5] 


[750; +CX)] 


-1.16 ±0.04 


-23.83 ±0.03 


0.75 ±0.08 




This work 


CMB 


[-25 


-21] 


[5000; 20000] 


-0.73 ±0.02 


-23.17 ±0.02 


1.11 ±0.02 


3.94 ±0.02 


This work 


CMB 


[-25 


-17] 


[750; 20000] 


-0.80 ±0.01 


-23.22 ±0.01 


1.13 ±0.02 


4.16 ±0.02 


This work 


LG 


[-25 


-17] 


[750; 20000] 


-0.86 ±0.01 


-23.24 ±0.01 


1.13 ±0.02 


4.25 ±0.02 


This work 


LG 


[-25 


-21] 


[5000; 20000] 


-0.76 ±0.02 


-23.18 ±0.01 


1.14 ±0.02 


4.02 ±0.02 


|6| > 10, X < 11.5 


CMB 


[-25 


-17] 


[300; 20000] 


-0.94 ±0.02 


-23.28 ±0.01 






1/Vmax fit 


CMB 


[-25 


-21] 


[750; 20000] 


-1.03 ±0.02 


-23.43 ±0.01 


0.85 ±0.06 


4.22 ±0.11 



from low surface brightness galaxies and the low surface 
brightness regions of catalogued galaxies is of little con- 
cern to us. 

(iv) The fitting method itself may also make a differ- 
ence. Our default parametric fit is pinned to the mag- 
nitude range where the formal Poisson errors are small- 
est, namely [—25, —22]. However, we have seen that sys- 
tematic effects can be important. As an alternative, we 
have taken the LF given by the l/Knax method, added 
in quadrature the statistical error bars and the fluctua- 
tions from the different sub catalogues, and fitted these 
data with a Schechter LF. As indicated in Table [21 we 
have obtained a steeper faint end slope and a brighte r 



Schechter LF, is 



M* , w h ich are in better a gree ment with Kochanek et al.l 
(l200lh . lCole et al.1 (l200lh and [Huchra et al.1 (l2005l ). but 
still discrepant with I Jones et al.1 (|2006l ) . 



We conclude that, given all of these systematics, our 
LFs are reasonably consistent with those that have been 
found previously. One aspect which can be improved is 
peculiar velocity corrections, but we postpone a fully self- 
consistent treatment of peculiar velocities and the LF 
determination to a future paper. 

We confirm that the bright end part of the LF does 
not seem to follow a Schecht er LF, as already seen by 
the GdfGRS (I Jon es et al]l2006l ). This effect is clearly seen 
in the SDSS, 2MRS and GdfGRS subsamples separately. 
The deviation becomes significant at Mr ^ —25, or 
two magnitudes brighter than M*, and is presumably 
due to brightest cluster gal axies, which have typical K- 
band magnitudes of ^ — 26 (|Lin Mohrll2004l ) and have 
long been known t o deviate from the extrapolat ion of a 
Schechter function iTremaine Richstond (| 19771 ). 

We may check the consistency of this LF with the 
number of galaxies in the 2M++ catalogue. We predict 
that the total number of galaxies of redshifts between 
the distances rmin(2;min) and rmax(^max), assuming the 



N 



/ $(M)dM / d^r/(r,f,M^in,M„ 



(23) 



We plot this function as a solid green line in Fig. (5] We 
also show the predicted number of galaxies in each bin 
of the grey histogram by a solid blue line. We see that 
the prediction in each redshift bin agrees well with the 
observed number of galaxies, but the total is off by ^2%. 
The difference comes both from the low luminosity part 
of the luminosity which is not adjusted because of our 
cut at cz > 5,000 km s~^ and the high luminosity part 
for which objects are not following a Schechter LF, as in 
Fig. 13 

In Table [21 we also give the mean luminosity density 
L as derived from Eq. ([221 • ^ is a lot less sensitive than 
n to the faint end of the luminosity function. As before, 
the errors are dominated by systematics due to the differ- 
ent corrections from peculiar velocities and the adequacy 
of the Schechter function to fit the observed luminosity 
function. Taking the average and computing the disper- 
sion in values for L for the four tests indicated in Table [21 
yields L = (4.09 ± 0.12)10^ L©. 



2.7 Weights 

Using the LF, we may now compute the appropriate 
weights to give to observed galaxies to account for incom- 
pleteness of the redshift catalogue. Our long-term goal 
is to reconstruct the dark matter density, under the as- 
sumption that galaxies trace the dark matter. There are 
several ways to link the galaxy density to the dark matter 
density: assuming that there is a linear relation between 
the two fields, one might consider number-weighting, 
in which the DM density is assumed to be related to 
the number-density of galaxies, or luminosity-weighting, 
which can serve as proxy for stellar mass, and so may 
be a better tracer of DM density. We will consider both 
of these schemes here. More complicated relationships. 



© 0000 RAS, MNRAS 000, 000-000 



10 G. Lavaux & M. J. Hudson 



for ex ample based on a halo model (jMarinoni HudsonI 
will be considered in a future paper. 
We compute number- weighting based on the fraction 
of observed galaxies: 



observed (j^) 
-^average 
/(r,r,Mmin,Mmax) 

/^^^^ ^o(M) dM 



(24) 



The weight applied to each galaxy is then 1 / /observed ('^)- 
This procedure is coinmon and has been used pre viously 
(e.g. iDavis k Hucto]ll982l : IPike Hudso^l2005l ) 

We follow a similar procedure for correcting the local 
luminosity density of galaxies by estimating how much 
light we are missing at the distance of each galaxy located 
at position r. The fraction 



/observed (■'^) 
1 



-f^observed (l*) 



^f(^)^M^L^M^.^ (^/, ^10)) (25) 

of luminosity with ^observed is the mean luminosity ex- 
pected to be observed in a small volume at position r, 
r = r/\r\ and Laverage the mean luminosity emitted by 
galaxies in the Universe. The value of L 

average IS 



J M 



L(M)$o(M) dM. 



(26) 



The weight to apply to each intrinsic luminosity of a 
galaxy is then 1 / /observed ('^)- This procedure has already 
also been used with success w i th ob s ervati o n and mock 
catalogues fe.e:. Lavaux et al. 1 I2OO8I . I2OIQI : iDavis et al.1 

[mil). 

For our choice of absolute magnitudes, Mmin = 
-25 and M^ax -20, the 2M++ is volume lim- 
ited up to rmin ^ 20h~^ Mpc, and extends up to 
T^max =S00h~^ Mpc. We find that at a distance of ^ 
150h~^ Mpc, the galaxy number weights are typically 
between 10 and 400, depending on whether the region is 
limited to K2M++ < 12.5 or K2M++ < 11-5, respectively. 
Similarly, the luminosity weights range between ^2 and 
~40. So weighing by luminosity has the advantage that 
it is less noisy at large distances. 



3 TREATING THE ZONE OF AVOIDANCE 

The "Zone of Avoidance" (ZoA) is the region of the 
Galactic plane where observations of galaxies are difficult 
due to the extinction by Galactic dust and stellar confu- 
sion. We show in Fig. [51 the number of galaxies in 2M++ 
with i^2M++ < 11-5, in bins of sin(6), and corrected for 
incompleteness effects. We see that the distribution is 
close to flat as a function of Galactic latitude, except for 
a hole contained between the latitudes —10° < b < 10°. 
We define the ZoA in 2M++ as this band for galactic 
longitudes —30 < / < 30, but reduce it to 5° outside this 
range. In addition, we impose the constraint that the ab- 
sorption not to exceed Ak = 0.25 in regions devoid of 
galaxies. 




0.0 

sin(6) 



Figure 8. The effect of the ZoA on 2M++. The weighed 
number density of galaxies in each bin of sin(6) is shown by 
the thin solid histogram. The dashed line shows the number 
density of galaxies once ZoA is filled with cloned galaxies. The 
two thick vertical lines correspond to 6 = ±10°. Here we used 
only galaxies for which cz < 15, 000 km s~^. 



In order to reconstruct the density field over the full 
sky, it is clearly necessary to fill the ZoA. One option 
is to fill it with mock galaxies so that their density of 
these objects matches the mean density outside the ZoA. 
This option, would however, fail to interpolate large-scale 
structure observed above and below the ZoA. The option 
adopted here, following iLvnden- Bell et all (|l989l ). is to 
"clone" galaxies immediately above and below the ZoA. 
The procedure of creating a galaxy clone at a latitude be 
of a galaxy at latitude b is simply to shift the latitude 



sm(bc) = sm(b 



- sm{b) 



(27) 



where 



if\b\ > 5°and|/ 
if 161 > 10° and I 



> 30° 
I < 30° 



(28) 



sign(6) X 5° 
sign(6) X 10° 

We refer to these as "ZoA-clones" . Fig. [8] shows the dis- 
tribution of galaxies as a function of Galactic latitude 
before and after cloning. After cloning, the distribution 
shows no dependence on latitude. 



4 GROUPING GALAXIES 

We use redshifts to estimate galaxy distances, but in the 
presence of peculiar velocites this relationship is not per- 
fect. In addition to the so-called ' Kaiser (1987) effect" 
which aflFects very large scales, there is also a contami- 
nation by the "Finger-of-god eflFect" due to the velocity 
dispersion of galaxies in clusters of galaxies. This causes 
a significant amount of noise on the redshift-estimated 
distance. One way to deal with the problem is to group 
galaxies, which by simple averaging improves the dis- 
tances estimated from the redshifts. The grouping in- 
formation is also interesting to study the statistics and 
properties of galaxy groups. 

In this Section, we describe the algorithm used to 
assign galaxies to groups and clusters. We use this infor- 
mation in the next sections for deriving a better density 
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10^ 



10^ 




^•^^0 5000 10000 15000 20000 

czcuB (km.s-^) 

Figure 9. Number of group members (richness) as a function 
of redshift. The richness is is not corrected for incompleteness. 



IQi^ 




10^ 



^^ 5000 10000 15000 20000 

c^cMB (km.s-^) 

Figure 10. Group luminosity as a function of redshift. The 
luminosities are not corrected for incompleteness. 

field fSection l5.ip and, in a future work, peculiar veloci- 
ties. Grouping also allows a better determination of the 
center of mass of superclusters (Section [521) their in- 
fall pattern (Section l5.3p . As a byproduct of 2M++, we 
provide the a catalog of groups and their properties in 
Appendix [BJ 

4.1 Grouping Algorithm 

To assign galaxies to groups we use the standard per- 
colation, or "Friends-of -friends" algorithm developed by 
iHuchra GeileJ (|l982l ). The algorithm is designed to 
identify cone-like structures in redshift space. Two galax- 
ies are considered to be part of the same group if: 

- their estimated angular distance separation is less 
than Dsep, 

- their apparent total velocity separation separation is 
less than Vq. 

Vq is kept fixed for the whole volume of the catalogue. 
Dsep is adapted such that the detected structures are 
always significant compared to the apparent local number 



Group name 




^sep 




b 


Virgo 
Fornax 


-oo 2,500 km s^^ 
-oo 1,600 km s-i 


10° 

8° 


279 
240 


74 
-50 



Table 3. Parameters for manual grouping of galaxies. All 
galaxies which are in the direction (/, b) and within the red- 
shifts [^^min; 2;max] with a maximum angular separation to (I, b) 
equal to ^sep, are considered part of the group indicated in the 
first column. 



density of galaxies, by explicitly accounting for selection 
effects. The constraint of a constant local overdensity at 
a redshift distance z leads to 

^sep = Do t-'-' ' , (29) 

with Lmin{z) the minimum absolute luminosity observ- 
able at redshift ^{L) the galaxy LF, zf the fiducial 
redshift, Dq the selection angular distance at fiducial red- 
shift. The parameter Do is linked to the sought overden- 
sity ^overdensity for group dctcctiou by the relation 

^overdensity = — ^^0 / ^L) dL ] - 1. (30) 

\ ^ JL^iAzp) J 

This density is computed at the fiducial redshift distance 
Zf. 

We have chosen the following parameters for defining 
our groups: Vf — czf — 1,000 km s~^ and ^overdensity = 
80. These parameter s have been used in previous studies 
([Ramella et al.|[l989l ). With the LF, for our choice of fidu- 
cial parameters, we compute that the transverse linking 
length is Dq = 0.45/i"^ Mpc. We count 4 002 groups with 
three or more members within 2M++, for redshift dis- 
tance less than 20,000 km s~^. We do not group galaxies 
farther than 20,000 km s""*^ where the catalogue becomes 
sparse. For the very nearby Virgo and Fornax clusters, 
the FoF algorithm fails and so we manually assign galax- 
ies to nearby clusters according to the parameters given 
m 

table El 



4.2 Results 

In Figure [9l we plot the richness of detected groups as a 
function of redshift. In Figure [TOl we have plotted the 
total luminosities of the same groups. Finally, in Fig- 
ure llli we give the velocity dispersion of the galaxies 
within these groups. As expected the mean velocity dis- 
persion does not vary significantly with distance, and has 
a mean value of ^ 95 km s~^. The richness is approxi- 
mately constant up to ~ 150/i~^ Mpc, which is a design 
feature of the group finder. The minimal luminosity of the 
groups increases with distance, as we are losing the fainter 
objects at larger distances because the 2M++ catalogue 
is limited in apparent magnitude. The catalogue of group 
properties is given in Appendix [B] We have checked that 
the parameters of the fitted Schechter LF do not change 
significantly after the grouping of the galaxies. 
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Shapley super cluster 



"^■^^H^rc^les cluster , 



Hydra-CentaurusV^ 
supercluster 



Norma cluster 




Virgo cluster 
Milky Way 



Perseus-Pisces 
supercluster 



Figure 13. The 2M++ galaxy distribution and density field in three dimensions. The cube frame is in Galactic coordinates. The 
Galactic plane cuts orthogonally through the middle of the back vertical red arrow. The length of a side of the cube is 200h~^ Mpc 
and is centred on Milky Way. We highlight the iso-surface of number fluctuation, smoothed with a Gaussian kernel of radius 
1,000km s"-*^. Si, = 2 with a shiny dark red surface. The position of some major structures in the Local Universe are indicated 
by labelled arrows. We do not show isosurfaces beyond a distance of 150h~^ Mpc, so Horologium- Reticulum is, for example, not 
present. 



10^ 




^°"o 5000 10000 15000 20000 

c^cMB {km.s-1) 

Figure 11. Group velocity dispersions as a function of red- 
shift. The thick solid line indicates the trend of the evolution 
of the average velocity dispersion with redshift. The scale of 
the variation is ~ lAh~^ Gpc, far larger than the depth of 
the 2M++ catalogue. 



5 DENSITY FIELD 

In this Section, we consider some properties of the peaks 
in the three-dimensional density field obtained from the 
distribution of galaxies in the 2M++ galaxy redshift cat- 
alogue. We assume that the number density and luminos- 
ity density of galaxies follow a Poisson distribution. As 
such, the mean smoothed density contrast p(x) given the 
galaxy weights Wi is 

TV 1 

-'^galaxies 

p(s) = - y V^(s-s.)^z (31) 
and the standard deviation 

-^galaxies 

4i^) = ^ E W{s-s,fwl (32) 

with s the coordinate in redshift space, VK(x) the smooth- 
ing kernel considered. To compute the position of the 
peaks in this density field, we use an iterative spherical 
overdensity algorithm: 

(i) we initialize the algorithm with an approximation 
Xc of the expected position of the cluster; 
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SGX (km.s-i) 

Figure 12. The 2MH — h number density field in Supergalactic 
plane. The density field is smoothed with a Gaussian kernel of 
1,000 km s"-*^ radius. Colour contours show the overdensity in 
units of the mean density and are separated by 0.5. The mean 
density is highlighted by a white contour. 

(ii) we compute the bary center x^^-^ of the set of 
galaxies contained in a sphere centred on and with 
radius Rn] 

(iii) we iterate (ii) until convergence, setting Rn-\-i = 

(iv) we reduce Rn+i = O.SORn- If Rn+i > lh~^ Mpc, 
then we go back to step (ii), in the other case we termi- 
nate the algorithm. 

We define the position of the structure as the one given 
by the last step in the above algorithm. This position is 
used in the following sections to compute mean densities 
and infall velocities on clusters. 



5 . 1 Cosmography 

Fig.ll2lshows the galaxy number density field of our cata- 
logue in the Supergalactic Plane, smoothed to 10h~^ Mpc 
with a Gaussian kernel. The Shapley concentration 
(SC) in the upper-left corner, near (SGX, SGY) ^ 
(—10000, 7000) km s~^, is particularly prominent and is 
the largest density fluctuation in the 2M++ catalogue. 
The Shapley region is covered by the 6dF portion of 
the survey which extends to a depth K2M++ < 12.5. 
Shapley is thus correctly sampled and is not a result of 
overcorrection of data limited to X2M++ < 11.5. When 
smoothed with a Gaussian kernel of 10h~^ Mpc radius, 
the Shapley concentration peaks at (/, 6) = (312, 30) and 
d = 152/i"^ Mpc with a density 1 + = 8.83 di 0.46, in 
galaxy number density contrast, and 1 + & = 9.51±0.54 
in terms of luminosity density contrast. 

The second most important structure in the Super- 
galactic plane of the 2M++ catalogue is the Perseus- 
Pisces (PP) supercluster. It is clearly seen in the 
Supergalactic plane in Fig. [12] at (SGX, SGY) 
(5000, -1000) km s"^ Its highest redshift space density. 



smoothed with Gaussian kernel of 10h~^ Mpc radius, is 
about 1 + = 4.46 db 0.18 in terms of number density 
contrast, 1 + = 4.47 ± 0.20 in terms of luminosity den- 
sity contrast. The position of the peak corresponds to 
the Perseus cluster at (/, 6) = (150, —13) , which is quite 
near the ZoA, and a distance of 52 Mpc. It is quite 
possible that the filling of the ZoA by galaxies cloned 
from the Perseus itself amplifies the overdensity of this 
supercluster. 

The extended overdense structure in the central part 
of the Supergalactic plane, at about (SGX, SGY) ^ 
(-5000,0) km s"\ is the Hydra- Cent aurus- Virgo (HC) 
supercluster. At 10h~^ Mpc smoothing scale, the high- 
est peak, located at (/,6) = (302,21), d = 38h~^ Mpc, 
coincides with the Centaurus cluster and has a height 
of 1 + = 3.02 ± 0.08 in number density contrast and 
1 + = 3.40 ±0.14 is luminosity density contrast. 

Finally, Fig.[T3]is a three-dimensional representation 
of the catalogue in Galactic coordinates, which means 
the Galactic plane goes through the middle of the ver- 
tical sides of the box, near the Norma cluster. We plot 
the 2M++ galaxies as points. Strong overdensities are 
highlighted by a transparent dark-red iso-surface of den- 
sity fluctuation of luminosity 6^ = 2. This density has 
been smoothed at 10h~^ Mpc with a Gaussian kernel 
from the corrected number distribution. The Shapley su- 
percluster is located at the top-left corner of the cube. 
A number of overdensities in the right part of the cube 
arise from the high weights, as this region has a depth of 
only i^2M++ < 11.5. 



5.2 Supercluster masses 

We show in Fig. [TJ] the mean overdensity and ex- 
cess mass within a sphere of 50 Mpc for four 
important superclusters in the 2M++ catalogue: the 
Shapley concentration, the Perseus-Pisces supercluster, 
the Horologium-Reticulum (HR) supercluster centred at 
(/,6) = (265,-51) and a distance of 193 Mpc and 
the Hydra- Centaurus supercluster. The profiles are cen- 
tred on the position where the density peaks for each 
supercluster. 

For the four superclusters, we note that the pro- 
files obtained through number weighing and luminosity 
weighing are nearly equivalent. The bumps in the mean 
density, shown in the left panels, are reproduced in both 
weighing schemes. This is particularly striking for the PP 
supercluster, even for scales as small as 10h~^ Mpc. In all 
cases, the luminosity weighted contrast is slightly lower 
than the number weighted density contrast. In the follow- 
ing discussion, we adopt the luminosity- weighted number 
contrast. 

The excess masses of all superclusters converge at 
radii of ^ 50 Mpc. While Shapley is the most mas- 
sive supercluster, we find that HR is very similar when 
measured on scales of 50 Mpc. Both have masses 
close to 10^^ h~^MQ. The PP and HC superclusters are 
less massive, but, being considerably closer, these have 
more impact in the motion of the LG and nearby galax- 
ies, as we discuss below. 

Our estimate of Shapley's mass and density contrast 
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Figure 14. The cumulative average density profile and the excess mass as a function of radius from four major superclusters 
in the 2MH — h redshift catalogue. In the two panels, we both show the profiles computed using the number weighed (solid lines) 
and the luminosity weighed (dashed lines) scheme. In the left panel, the horizontal black dashed line corresponds to the mean 
density. In the right panel, the dotted lines indicates the mass of a sphere of the given radius at the mean density. Note that 
since we are plotting excess mass, to obtain the total mass one must add this value. The black, dark grey, light grey and red lines 
correspond respectively to the Shapley concentration, the Horologium-Reticulum supercluster, the Perseus-Pisces supercluster and 
the Hydra-Centaurus supercluster. The error bars are estimated assuming that galaxies follow Poisson distribution for sampling 
the matter density field, as given by Eq. (|32|) . 



is similar to that of IProust et al.l (|2006l ) who measured 
a density contrast 5n — 5.4 di 0.2 in a truncated cone of 
225 square degrees between 90 and 180 Mpc with 
a volume equivalent to a sphere of effective radius 30.3 
Mpc. In a sphere of this radius centred on Shapley, 
we f ind a luminosity densit y contrast of 5k,l = 4.1 ±0.15. 

iMufioz Loebl (|2008l ) calculated the mass of SC 
based on the overdensity of rich clusters and obtained 
a mass 3.3 ± 0.3 x 10^^ h~'^MQ within a sphere of 
35 Mpc. On the same scale, we obtain a mass of 
4.87 ± 0.18 X 10^^ /i"^M0, assuming = 0.3 and 
bK,L = 1 for K-hand l umino sity. Using similar argu- 
ments, [Sheth^maferi^ (|201lh quote a mass of 1.8 x 10^^ 
/i~^M0 within a slightly smaller radius of 31 Mpc. 
On the same scale we find 4.00 ± 0.17 x 10^^ /i^^M©. 
These values could be brought into rough agreement 
if luminosity- weighted 2MASS galaxies are strongly bi- 
ased, with b = 2-3. Such a strong biasing would, 
however, conflict wi th the measurement b K,n — 1-05 ± 
0.10(Qm/0.3)°-^^ bv IPike HudsonI (l2005l ) but may be 
marginally consistent with the lowe r value bK,n = 1.56 =b 
0.16(Qm/0.3)°-^^ found recently bv lPavis et aP (|201lh . 

A further caveat is that our density estimates are 
in redshift -space, and s o are enhanced by a factor up to 
ps — 1.2 (iKaisej 1987l lFI c ompared to the estimate s of 
iMufioz Loebl (|2008l ) and ISheth Diaferiol (|201lh . In 
a future paper, we will reconstruct the density field in 
real-space and calibrate the biasing factor directly using 
peculiar velocity data, so a detailed comparison of over- 
den sities awaits fu t ure w ork . 

[Hudson et al.l (|2004l l studied the overdensity of the 
SC as traced by IRAS-selected galaxies. Within a 50 



10 



Shapley 

Horologium 

Perseus 

Hydra-Centaurus 




Figure 15. The infall velocities as a function of distance 
for four major superclusters in the 2M++ redshift catalogue. 
Curves and error bars are as in Fig. fl4l 



V1 + 2//3 + /2/5, with / = 0.5. 



Mpc-radius sphere they found that the overdensity 
of IRAS-selected galaxies is only 0.2. Here we find that 
the overdensity of 2MASS-selected galaxies on the same 
scale is ^ 1. Clearly, the relationship between IRAS and 
2MASS-selected galaxies is not well-described by a rela- 
tive linear bias, since a value of ^ 5 would be required in 
the Shapley supercluster, whereas the field requires a rel- 
ative b ias between 2MASS- and IRAS-selected galaxies 
of - 1 (IPike Hudsonll2005l ). 
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Supercluster 


Sphere centre 

7 


b 


1 + 


Mass 




Mpc 


o 


o 




10^6 H-^Mq 


Shapley 


152 


312 


30 


2.05 ±0.05 


8.9 ±0.2 


Horologium- Reticulum 


193 


265 


-51 


2.01 ±0.10 


8.7 ±0.5 


Perseus-Pisces 


52 


150 


-13 


1.41 ±0.03 


6.1 ±0.1 


Hydra-Centaurus 


38 


302 


21 


1.43 ±0.03 


6.2 ±0.1 



Table 4. Luminosity density contrast and estimated masses (assuming Qm = 0.3 and bK,L = 1) of the four superclusters from 
the distribution of galaxy light within a sphere of radius 50 Mpc. 



5.3 Supercluster infall 

We now discuss the impact these structures have on large- 
scale flows in the nearby Universe. We have estimated the 
infall velocity onto each of these structures using linear 
theory: 

^^infaii = ^13H6{R)R, (33) 

with H the Hubble constant, /3 = f/b where / is the 
linear density perturbation growth rate and 6 is a biasing 
parameter, and S{R) the mean density inside a sphere of 
radius R and centred on the supercluster. For a ACDM 
cosmology, / ^ (|Bouchet et al.lll995|V We use (3 = 
0.5 (as determined for 2MASS galaxies by Pike HudsonI 
[2005i) whenever we need to estimate a velocity. This value 
corresponds to Qm — 0.30 and Qa = 1 — Qm with b = 1. 

Fig. [15] shows the infall velocity profiles of the four 
superclusters. Although we plot the linear theory infall 
down to small radii (i?<10/i~^ Mpc), we note that lin- 
ear theory does not apply in these regions and focus the 
discussion on distances 10h~^ Mpc. The infall veloc- 
ities at 10h~^ Mpc are all at least 2,000 km s~^, with 
Shapley having the highest infall at nearly 4,000 km s~^. 
At 50h~^ Mpc, the Shapley and HR superclusters have 
an infall of '^SOO km s~^. The average overdensity of the 
Shapley concentration within a sphere of 50h~^ Mpc is 
1 + = 2.05 ± 0.05. Neglecting structures beyond 50 

Mpc, linear theory implies that the supercluster is 
responsible for attracting the LG with a peculiar veloc- 
ity of 90 ± 10km s~^. This motion represents ^ 15% of 
the total velocity of the LG with respect to the CMB 
rest frame. Although the excess mass of HR is similar to 
Shapley, its effect on the LG's motion is less than that 
of Shapley due to its greater distance: we estimate 60 
km s~^. Added vectorially, the net peculiar velocity from 
these two superclusters is approximately 110 km s~^ to- 
wards (/,6) = (297, —1). This direction is within the er- 
ro rs of the direct i on of the 407 km s~^ bulk flow found 
bv lWatkins et al.1 (|2009l \ but is lower in amplitude. 

Because they are closer to the LG, the HC and PP 
superclusters have a greater impact. Approximating HC 
as a sphere, the infall at the LG's distance of 38 Mpc 
is 588 ± 26km s~^. Whereas PP is denser, its greater 
distance of 52 Mpc puts it on the losing side of the 
gravitational tug-of-war with HC: the infall of the LG 
towards PP is only 313±24 km s~^. Thus the net motion 
is towards HC. 

Note that it is likely that underdense regions also 



contribute a push. iKocevski Ebelind (|2006bl ) have 
noted the deficit of rich clusters in the Northern sky, par- 
ticularly in the distance range 130 to 180 Mpc. Thus 
a full analysis of peculiar velocities requires integration 
over the entire density field, a topic we defer to a later 
paper. 



6 SUMMARY 

We have compiled a new, nearly full- sky galaxy redshift 
catalogue, dubbed 2M++, based on the data from three 
redshift surveys: the 2MASS Redshift Survey {K <11.5), 
the Sloan Digital Sky Survey and the 6dF galaxy redshift 
survey. After having calculated corrected magnitudes and 
having calculated redshift completeness, we have deter- 
mined LFs and weights that allow us to determine the 
redshift density field to a depth of 200 Mpc. The 
most prominent structure within 200 Mpc is the 
Shapley Concentration: its luminosity density within a 
sphere of radius 50h~^ Mpc is 2.05 times the mean, and 
is thus responsible for approximately 90 km s~^ of the 
LG's motion with respect to the CMB rest frame. We 
have compared the density profile of four massive su- 
perclusters that are present in the 2M++ catalogue: the 
Shapley Concentration, the Perseus-Pisces supercluster 
and the Horologium-Reticulum supercluster and Hydra- 
Centaurus. The Shapley Concentration is clearly the 
most massive of the four, but HR is only slightly less 
massive. 

This new, deep full-sky catalogue will be used in fu- 
ture work to study the peculiar velocity of the LG and 
other nearby galaxies. Our hope is that the distribution 
of density in the 2M++ volume will account for the 
high- amplitude bulk motions on scales of 100h~^ Mpc 
(|Watkins et al.1 1 20091 : [Lavaux et al.ll201Ql ). 
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none 


03355460-8537067 


299.59 


-30.39 


11.90 


12702 


12736 


270 


3775 


1.0 


0.7 





1 








1 


none 


08423963-8430223 


297.59 


-24.47 


11.70 


12284 


12357 


270 


4039 


1.0 


1.0 





1 








1 


none 


08431792-8429053 


297.58 


-24.44 


12.35 


12284 


12357 


270 


4039 


1.0 


1.0 





1 








1 


none 


08424060-8427453 


297.55 


-24.44 


11.96 


12284 


12357 


270 


4039 


1.0 


1.0 





1 








1 


none 


ZOAOOOOOOO 


330.07 


-0.41 


10.21 


5238 


5387 
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0.7 


1 
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zoa 


ZOAOOOOOOl 


330.49 


-0.33 


11.31 


8841 


8989 
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0.6 
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zoa 


ZOA0000002 


330.01 


-1.18 


11.96 


11830 


11983 







1.0 


0.9 


1 
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zoa 


ZOA0000003 


331.29 


-1.10 


11.09 


3158 


3306 


70 




0.9 


0.8 


1 











1 


zoa 


07243410-8543223 


298.20 


-26.43 


11.51 


5301 


5361 


69 


4638 


1.0 


0.7 














1 


20096dF 


03403012-8540119 


299.56 


-30.29 


12.41 


12714 


12749 


66 


3775 


1.0 


0.6 














1 


20096dF 


07400785-8539307 


298.20 


-26.13 


11.09 


5184 


5246 


11 


4638 


1.0 


0.7 














1 


2008ApJ 


03355460-8537067 


299.59 


-30.39 


11.90 


12702 


12736 


270 


3775 


1.0 


0.7 





1 








1 


none 


07420104-8525161 


297.96 


-26.04 


10.31 


5150 


5213 


14 


4638 


1.0 


0.7 














1 


20096dF 


02090195-8520255 


301.12 


-31.51 


11.52 


12675 


12699 


11 


3776 


1.0 


0.9 
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20112MRS.... 



Table Al. The 2M++ catalogue - Col. (1): the name of the galaxy as given in the 2MASS-XSC database. Col. (2): Galactic longitude in degrees. Col. (3): Galactic latitude in degrees. 
Col. (4) : Apparen t mag nitud e in band Ks a^s d efined in Section [2?2] Col. (5): Heliocentric total apparent velocity. Col. (6): Total apparent velocity in CMB rest frame, using relation 
from lKogut et al.l (Il993l ) and lTullv eraP (l2008l V Col. (7): Total apparent velocity error (equal to zero if not measured). Col. (8): Unique group identifier obtained from the algorithm of 
Section m Col. (9): Redshift incompleteness at magnitude K2M++ ^ 11.5. Col. (10): Redshift incompleteness at magnitude K2M++ < 12.5. It may be empty in that case the catalogue 
is limited to X2M++ ^ 11-5 in the portion of the sky holding the galaxy. Col. (11): Flag to indicate whether this is is a fake galaxy to fill the zone of avoidance following the algorithm 
of Section O Col. (12): Flag to indicate if the redshift has been obtained by the cloning procedure of Section [231 Col. (13): Flag to indicate whether this galaxy lies in the exclusive 
region covered by the 2MRS target mask (2Mx6S region). Col. (14): Flag to indicate whether this galaxy lies in the non-exclusion region covered by the SDSS. Col. (15): Same as (14) 
but for the GdFGRS. Col. (16): Bibliographic code for the origin of the redshift information. The code is truncated in the above table but available in full in the electronic version of 
the catalog. 
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UJ 




1 


281.26 


73.47 


6.60 


65 


1223 


1556 


599 


1000 


182.41 


-13.06 


10.65 


4 


3 


-25 


97 


1001 


91.75 


51.01 


6.15 


4 


691 


752 


51 


1002 


137.70 


12.33 


5.39 


9 


1139 


1055 


118 


1003 


316.23 


-10.88 


10.42 


3 


-80 


17 


21 


1004 


123.60 


74.51 


3.80 


9 


-249 


-33 


156 


1005 


184.67 


83.05 


4.96 


24 


675 


956 


312 


1006 


144.05 


66.22 


3.59 


69 


845 


1049 


227 


1007 


171.51 


32.81 


8.58 


4 


506 


646 


46 


1008 


108.51 


58.06 


5.83 


6 


186 


302 


109 


1009 


319.07 


-12.21 


4.75 


4 


-146 


-65 


89 


1010 


144.70 


36.20 


3.55 


5 


92 


157 


82 


1011 


33.58 


14.01 


8.21 


4 


1865 


1780 


101 


1012 


41.38 


14.94 


8.00 


3 


2291 


2188 


12 


1013 


134.90 


32.70 


6.94 


4 


1330 


1348 


51 


1014 


103.67 


33.03 


8.10 


4 


1174 


1129 


63 


1015 


41.06 


12.68 


9.01 


3 


2741 


2626 


93 


1016 


45.04 


17.71 


8.24 


3 


2261 


2162 


60 


1017 


150.98 


5.93 


8.69 


4 


5096 


5028 


112 


1018 


147.84 


7.90 


8.10 


6 


4807 


4736 


237 


1019 


129.32 


8.92 


7.37 


5 


3276 


3146 


203 


1020 


103.21 


12.39 


7.94 


5 


2692 


2523 


131 


1021 


50.91 


6.89 


8.34 


4 


4831 


4659 


193 


1022 


69.13 


8.13 


8.96 


4 


4556 


4359 


81 


1023 


75.62 


6.03 


8.22 


4 


4711 


4497 


94 


1024 


269.12 


5.57 


9.91 


3 


5052 


5324 


75 


1025 


264.14 


7.22 


8.89 


4 


4685 


4965 


234 


1026 


275.43 


8.94 


9.48 


3 


4008 


4291 


88 


1027 


264.36 


8.40 


8.29 


7 


4796 


5081 


121 


1028 


295.46 


8.81 


8.36 


4 


4449 


4701 


113 



Table Bl. The 2M++ group catalogue - Col. (1): Group identifier in the catalogue. It corresponds to column 8 of Table IaT] Col. (2): Galactic longitude Col. (3): Galactic latitude 
Col. (4): Apparent magnitude in band Kq, as defined in Section [2?2] The magnitude is derived from the 2M++ galaxies. This is a magnitude uncorrected for incompletenes s effect. Col. 
(5): Richne ss, uncorrec t ed for incompleteness effect. Col. (6): Heliocentric total apparent velocity. Col. (7): Total apparent velocity in CMB rest frame, using relation from lKogut et all 
11993,) and lTullv erZI (120081) . Col. (8): Velocity dispersion in the group. 



